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ABSTRACT 

Emerging sonography techniques often imply increasing in the num- 
ber of transducer elements involved in the imaging process. Conse- 
quently, larger amounts of data must be acquired and processed by 
the beamformer. The significant growth in the amounts of data ef- 
fects both machinery size and power consumption. Within the clas- 
sical sampling framework, state of the art systems reduce processing 
rates by exploiting the bandpass bandwidth of the detected signals. 
It has been recently shown, that a much more significant sample- 
rate reduction may be obtained, by treating ultrasound signals within 
the Finite Rate of Innovation framework. These ideas follow the 
spirit of Xampling, which combines classic methods from sampling 
theory with recent developments in Compressed Sensing. Applying 
such low-rate sampling schemes to individual transducer elements, 
which detect energy reflected from biological tissues, is limited by 
the noisy nature of the signals. This often results in erroneous param- 
eter extraction, bringing forward the need to enhance the SNR of the 
low-rate samples. In our work, we manage to achieve such SNR en- 
hancement, by beamforming the sub-Nyquist samples obtained from 
multiple elements. We refer to this process as "compressed beam- 
forming". Applying it to cardiac ultrasound data, we successfully 
image macroscopic perturbations, while achieving a nearly eight- 
fold reduction in sample-rate, compared to standard techniques. 

Index Terms — Array Processing, Beamforming, Compressed 
Sensing (CS), Finite Rate of Innovation (FRI), Ultrasound, Xam- 
pling 

1. INTRODUCTION 

Diagnostic sonography allows visualization of body tissues, by 
radiating them with acoustic energy pulses, transmitted from a 
transducer array. Applying appropriate delays to the transmitting el- 
ements, in a process known as "beamforming", the interfering wave 
pattern forms a narrow beam, along which most energy propagates. 
As it propagates in the tissue, echoes are scattered by density and 
propagation- velocity perturbations [Tj|. These reflections are de- 
tected by the transducer array in a reception phase, where a second 
beamforming process is applied, aimed at localizing the scatter- 
ing structures while improving Signal to Noise Ratio (SNR) Q. 
The reception beamforming process involves averaging the detected 
signals, after aligning them with time-varying delays. Ultrasound 
systems perform this process digitally, requiring that the analog 
signals first be sampled. Confined to classic Nyquist-Shannon sam- 
pling theorem |3 |, these systems must sample the signals at twice 
their baseband bandwidth, in order to avoid aliasing. 

The goal of our work is to produce ultrasound images from sam- 
ples taken at rates far below these dictated by the classical frame- 
work. A first step in this direction was proposed in |4|, where the 
signals detected in each element were assumed to satisfy the FRI 



property 0. Applying the FRI framework, (4l|6) derive schemes, 
which allow reconstruction of ultrasound signals from samples taken 
at sub-Nyquist rates. These approaches follow the spirit of Xam- 
pling 1 7 1, a term which refers to the combination of classical sam- 
pling theory with recent developments in Compressed Sensing (CS). 
However, applying such schemes to data reflected from biological 
tissues, the resulting signal parameters are often erroneous. This is 
mainly due to the noisy nature of multiple reflections, also known as 
"speckles" Q, which are scattered by microscopic perturbations in 
the tissue. 

In order to cope with the noisy nature of the signals detected 
in individual elements, 1 8 1 first introduces a multi-sensor Xampling 
scheme, which enhances the SNR of the low-rate samples. This 
is achieved, by integrating traditional beamforming with FRI Xam- 
pling. The approach proposed in 1 8 1 is farther developed in (9), re- 
sulting in a more straightforward sampling mechanism. Both ap- 
proaches, referred to by the term "Compressed Beamforming", are 
derived and simulated for the special case of linear, B-mode scan. 

In the current work, we generalize the method introduced in 1 8 
9 1 to the case of polar scan. Additionally, we propose reconstructing 
the sampled signal by applying CS techniques, rather than spectral 
analysis methods, traditionally used in the FRI framework |5|. Fi- 
nally, we apply our approach to real cardiac ultrasound data, where 
we successfully image macroscopic perturbations, while achieving 
a nearly eight-fold reduction in sampling rate, compared to standard 
imaging techniques. 

This paper is organized as follows: in Section [2] we outline the 
principles of beamforming in ultrasound imaging. In Section [3] we 
discuss the FRI model in the context of ultrasound imaging, and its 
contribution to rate reduction when sampling signals at individual 
receivers. Our multiple-receiver Xampling scheme, which integrates 
beamforming into the Xampling process, is presented in Section [4] 
In Section[5]we introduce experimental results, obtained for cardiac 
ultrasound data. 

2. BEAMFORMING IN ULTRASOUND IMAGING 

We begin by outlining the beamforming process, lying at the heart 
of ultrasound imaging. The process is carried out both in trans- 
mission and reception. During transmission, modulated pulses are 
transmitted from all active array elements. Applying appropriate de- 
lays to the transmitting elements, the propagating energy may be 
directed along a narrow beam, in any desired direction, the result 
of a constructive and destructive interference pattern. The acoustic 
reciprocity theorem allows to perform a similar process during re- 
ception. In this case, we regard the perturbations, which reflect the 
transmitted energy, as point transmitters. Applying appropriate de- 
lays to the detected signals, the array may be focused to any coordi- 
nate, such that if, indeed, energy was scattered from that point, then 
the combined signals will undergo constructive interference. Know- 
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Fig. 1. Imaging setup. 



ing that most energy propagates along a narrow beam, and knowing 
its propagation velocity, we may associate time instances to radiated 
points in the medium, thus optimizing the reception focal point at 
these instances. The resulting process, in which time- varying delays 
are applied to the detected signals, prior to combining them for SNR 
enhancement, is referred to as "Dynamic Focusing". 

In order to derive an explicit expression for dynamic focusing, 
consider an array of M receivers, aligned along the x axis. We de- 
note by 5m the distance from the ra'th receiver to the origin. Assume 
a pulse of energy, transmitted at £ — from the origin, at direction 
0. The setup is depicted in Fig. [T] At time £, the pulse, propagat- 
ing at velocity c, crosses (x, z) — (c£cos#, ctsm6). A reflection 
originating from this coordinate will reach the ra'th element at time: 

r m (t) = t+ ^(tcosO) 2 + ( 7m - £sin<9) 2 , (1) 

with 7 m = 5 m / 1 c. Denoting by rao the index of a reference receiver, 
satisfying 5 mo = 0, it is readily seen that T mo (£) = 2t. Hence, in 
order to align the reflection detected in the ra'th receiver with the 
one detected in the reference receiver, we need to apply a delay to 
the detected signal, ip m (£), such that the resulting signal, (p m (£;#), 
satisfies (pm (2£; 6) — ip m (r m (£)). This implies a time-varying 
time delay, applied to cp m (£), resulting in the following distorted 
signal: 

(Pm(t;0) = (pm Q (t + v^ 2 - 47m£ sin + 4 7m ^ . (2) 
The aligned signals may be combined into the beamformed signal 

1 M 

*M) = ^Z>m(*;0), (3) 

m—l 

which exhibits enhanced SNR compared to {(p m (£; #)}^f =1 . 

Although defined over continuous time, ultrasound systems per- 
form the process formulated in 0>-(3| m me digital domain, requir- 
ing that the analog signals {(p m (£)} m=1 first be sampled. Confined 
to classic Nyquist-Shannon sampling theorem, these systems sam- 
ple the signals at twice their baseband bandwidth, in order to avoid 
aliasing. The detected signals typically occupy only a portion of 
their baseband bandwidth. Exploiting this fact, some state of the art 
systems manage to reduce the amount of samples transmitted from 
the front-end, by down- sampling the data, after demodulation and 
low-pass filtering. However, since such operations are carried out 
digitally, the preliminary sampling-rate remains unchanged. 

3. SAMPLE RATE REDUCTION USING THE FRI MODEL 

Our goal is to produce an ultrasound image from samples taken far 
below the Nyquist rate. A first step in this direction was proposed in 



|H , where the signals detected in each element were assumed to sat- 
isfy the FRI property. More specifically, they were approximated as 
L replicas of a known-shape pulse, h (£), caused by scattering of the 
transmitted pulse from L reflectors, located along the narrow trans- 
mission beam. Explicitly, the signal detected in the ra'th receiver 
was modeled as 

L 

Vrn (£) - ai,mh (t ~ £z,m), (4) 

1=1 

where £j jTn denotes the time, in which the reflection, originating 
from the /'th reflector, arrives at the ra'th receiver, and a^ m de- 
notes the reflection amplitude. Hence, the detected signal is com- 
pletely defined by 2L parameters, {£j, m , Q>i,m}f =1 - As shown in |5 |, 
spectral analysis techniques (e.g. annihilating filter [ 10|) may be 
applied, in order to extract the signal parameters, from its projec- 
tion onto a 2L-dimensional subspace, corresponding to its Fourier 
series coefficients. The sampling scheme is thus reduced to the 
problem of extracting K frequency samples of the detected signal, 
where K > 2L. Denote the temporal support of the detected signal, 
Lpm (£), by [0, T). Imaging to a nominal depth of 16cm, T is ap- 
proximately 210/isec. The signal's baseband bandwidth requires a 
nominal sampling rate of f s — 16Mhz, resulting in an overall num- 
ber of Tf s = 3360 real-valued samples. Still within classical sam- 
pling framework, assuming that the signal's passband bandwidth is 
only 4MHz, the data sampled at Nyquist rate may be down- sampled 
to 1680 real-valued samples. Thus, as long as the number of re- 
flected pulses satisfies L <^ 1680, substantial rate reduction may be 
achieved by an FRI Xampling scheme. 

The scheme proposed in (4j enables sub-Nyquist sampling of 
the signal detected in an individual array element. Sampling signals 
reflected from phantom targets, it managed to achieve good signal 
reconstruction. However, when applied to data reflected from bio- 
logical tissues, the resulting signal parameters are often erroneous. 
This is mainly the result of multiple reflections, scattered from mi- 
croscopic perturbations, manifested as noise. In extreme scenarios, 
where the noise masks the FRI component, the extracted parameters 
will be meaningless, such that any attempt to cope with errors in the 
parametric space will turn out useless. 

The novelty of our approach regards the assumption, that the 
FRI scheme applied to an individual element may be equally ap- 
plied to the beamformed signal defined in (3]), which obviously ex- 
hibits much better SNR. Indeed, it may be shown that, if the sig- 
nals {(fm (t)}m=i °^ e y me FRI model {4}, then $ (£; 6) is approx- 
imately of the form: 

L 

*(t;0) = 5>fc(*-*j)> (5) 

1 = 1 

and may thus be regarded as FRI. The parameters {ti, bi}f =1 in {5J 
depend on 0, yet we omit this dependency from our notation, for 
convenience. An obvious problem is that, as stated in Section [2] 
the signal <3> (£; 0) does not exist in the analog domain. In the next 
section, we derive our scheme, which manages to estimate its neces- 
sary low-rate samples from low-rate samples of filtered versions of 

4. COMPRESSED BEAMFORMING SCHEME 

The FRI Xampling scheme applies spectral analysis techniques, in 
order to estimate the 2L unknown signal parameters, from a subset 



of its consecutive Fourier coefficients. Denoting the cardinality of 
this subset by K, perfect reconstruction requires that K > 2L. 

Denote by k the set of indices {kj }f =1 , corresponding to the se- 
lected Fourier coefficients. Having obtained these coefficients, other 
mechanisms may be applied, in order to reconstruct the signal. In our 
work, for instance, we use CS methodology. Since each of the sig- 
nals {ifm (t)}™ =1 is supported on [0, T), then, by $2$-^ $ (£; 0) 
may only be constructed on a support contained in [0, T). Consid- 
ering its Fourier series, calculated with respect to this interval, let us 
explicitly write the fc/th Fourier coefficient, Cj. We will addition- 
ally assume, that the time delays, {£z}£ =1 , are quantized with A s 
quantization step, such that U = qiA s , q\ £ Z: 

c i = h r^hhit-^e-^^dt 

1 J ° 1=1 

= ^ V h f T h(t- qi A 3 )e- j ^ k ^dt 

1 1=1 J o 

= ±H(p j )± bl e-^^. (6) 

^ ' 1 = 1 

Here H (u) denotes the CTFT of h (£). Let N be the ratio [T/A s J . 
Then ([6} may be brought into the following matrix form: 



Here I[ a ,b) (t) denotes an indicator function having the value 1 for 
a < t < b and otherwise, and T m (0) is the maximal time value 
which is mapped, under the temporal distortion defined in (2), to the 
support of the beamformed signal. 

The process defined in |T0])-|TT]) can be translated into a 
multi-channel Xampling scheme, in which each signal, ip m (£), 
is multiplied by a bank of kernels {gj, m (£; @)}f = i defined by 
{TT}, and integrated over [0, T). This results in a vector c m = 
[ ci, m c 2 ,m ••• c K ,m } T . The vectors {c m }^ =1 are then 

averaged in c = [ a C2 ... ck ] T , which has the desired 
improved SNR property, and provides a basis for extracting the 
2L parameters which define (t; 6). Nevertheless, the complexity 
of the analog kernels, together with their dependency on 0, makes 
hardware implementation of such a scheme extremely complicated. 
We thereby take an additional step, approximating the samples 
{cj, m }f =1 from low-rate samples of <p m (£), taken in a much more 
straightforward manner. 

First, we substitute ip m (t) of |To] > by its Fourier series, calcu- 
lated with respect to [0, T). Denoting the n'th Fourier coefficient by 



[nj, we get: 



-HAx, 



(7) 



1 J ° 

^ 0m [kj - Tl] Qj,m;8 H, 



(12) 



where H is the K x K diagonal matrix with H (^-kj) as its j'th 
element, and x is a length N vector, whose j'th element equals bi 
for j — qu and otherwise. Finally, A is a K x N matrix, formed 
by taking the set k of rows from a N x N FFT matrix. Assuming 
that L <C N, this is a CS problem, aimed at estimating the L- sparse 
vector x from K samples, taken with a measurement matrix ^HA. 
A common technique for recovering x is the Orthogonal Matching 
Pursuit (OMP) algorithm, which we use in our work. 

We now return to the problem of obtaining the set k of $ (t; 6) 's 
Fourier coefficients. Let Cj denote the j'th Fourier coefficient of 
$(t;0). Then, 



° J ~ M ^ T J 

m=l J0 

where, from |2|: 



1 M 

'ipm (t;0)dt = — ^ 



(8) 



(9) 

Applying a change of variable to the integral,we can write ([9} as 



9j,m (t;0)ipm (t) dt, 



(10) 



with 



g jim (t]0) =q j ,m(t;6)e- i ^ k i\ 

qj,m(t]0) =I[\<y m \,T m (d)) (t) ^1 + 



7 m cos 2 
(t-jmSmOfJ " (11) 



. . 2ty . 7 m - t sin 

exp <; i — kj- —7m 



T "t- 7 m sin/: 



where Qj, m -,e [n] are the Fourier series coefficients of qj,m(t]6), 
also defined on [0, T). Replacing the infinite summation of |T2] > by 
its finite approximation: 



N 2 



Cj, m — ^ 0m [kj - n]Qj,m;0 N, 
n=Nx 



(13) 



it can be shown that, for any e > 0, and for any selection (j, ra; 0), 
there exist Ni (e, kj,m; 0) and A^2 (e, kj,m; 0) such that |cj, m — 
Cj,m\ < e- Therefore, with the finite sum |T3] >, we may approximate 
Cj, m to any accuracy. Furthermore, setting an upper bound on the 
energy of (p m (£), Ni and N2 may be chosen off-line, subject to the 
decay properties of the sequence {Qj,m-9 [n]} < ^L_ oo . 

With this result we can, for each and receiver index m, find 
the minimal set of ip m (t)'s Fourier series coefficients, from which 
we can simultaneously approximate the set {cj, m }f =1 . Denoting the 
set of indices corresponding to these coefficients by « m , arrange the 
set {4>m [k]} keKrn m a l en g m Km vector <I> m . Using |T3] >, we may 
define a linear transformation A m (0), such that: 



Cm A m (0) $ r 



(14) 



where c m is the length K vector containing Cj, m as its j'th element. 
Instead of extracting only K samples from each element using the 
complicated analog kernels, we now extract a larger number of sam- 
ples, K m . As we show in our simulation, good approximation is 
obtained with just a small sampling overhead. Our final Xampling 
scheme is illustrated in Fig. [2] Note the simple mechanism used 
for obtaining the Fourier coefficients in each individual element: a 
linear transformation, V m , is applied to point- wise samples of the 
signal, taken at a sub-Nyquist rate, after filtering it with an appropri- 
ate kernel, s m (— £), such as the Sum of Sines (4j. 
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(a) (b) (c) 

Fig. 3. Cardiac images generated by Xampling and using traditional methods, (a) standard beamforming applied to data sampled at Nyquist 
rate, (b) applying the non- approximated Xampling scheme, defined in (T0])-(pT|). (c) applying the final Xampling scheme of Fig. [2] 




Fig. 2. Xampling scheme utilizing Fourier samples of detected sig- 
nals. 



5. SIMULATION ON CARDIAC ULTRASOUND DATA 

We examine the result of applying our scheme of Fig. [2] to raw RF 
data, acquired and stored for cardiac images of a healthy consent- 
ing volunteer. The acquisition was performed using a breadboard 
ultrasonic scanner of 64 acquisition channels. The transducer em- 
ployed was a 64-element phased array probe, with 2.5MHz central 
frequency, operating in second harmonic imaging mode: 3 half cycle 
pulses are transmitted at 1.7MHz, resulting in a signal characterized 
by a rather narrow bandpass bandwidth, centered at 1.7MHz. The 
corresponding second harmonic signal, centered at 3.4MHz, is then 
acquired. The signal detected in each acquisition channel is ampli- 
fied and digitized at a sampling-rate of 50MHz. Data from all acqui- 
sition channels were acquired along 120 beams, forming a 60° sec- 
tor, where imaging to a depth of z = 16cm, we have T = 207/xsec. 
The results are illustrated in Fig. [3] 

The first image (a) was generated using the standard technique, 
applying beamforming to data first sampled at the Nyquist rate, and 
then down-sampled, exploiting its limited essential bandwidth. For 
a single image line, sampling at 50MHz, we acquire 10389 real- 
valued samples from each element, which are then down- sampled, 
to 1662 real-valued samples, used for beamforming. The result- 
ing image is used as reference, which we aim at reproducing with 
our Xampling scheme. We begin by applying our scheme without 
approximation, simulating the modulation with the analog kernels 
defined in Assuming L = 25 reflectors, and using two-fold 
oversampling, k comprises K — 100 indices. Since each sample 
is complex, we get an eight-fold reduction in sample-rate. The 
resulting image (b) well depicts strong perturbations observed in 
(a). Isolated reflectors at the proximity of the array (z « 6cm) 
remain well in focus. Applying our approximated scheme, for every 
kj G «, 1 < m < M and 0, we choose Ni and N2 of |T3] > such 



that \Qi,m-,e (^n)| 2 « 0.95 £~ -co |Qw (f n)\ 2 . 

At the current stage, this choice, and consequently the choice of 
{Km}m=i an d {A m }^ =1 , are performed off-line, based on the 
imaging setup. This results in an average number of 116 complex 
samples per receiving element, required for constructing each image 
line (with at most 133 samples required for some combinations of 
and element indices). We arrive at a seven-fold reduction in sample 
rate, and the resulting image (c) appears very similar to (b). 
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